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Abstract 

We derive the cross section of neutrino-nucleon scatterings in supernova magnetic fields, includ- 
ing weak-magnetism and recoil corrections. Since the weak interaction violates the parity, the 
scattering cross section asymmetrically depends on the directions of the neutrino momenta to the 
magnetic field; the origin of pulsar kicks may be explained by the mechanism. An asymmetric 
. neutrino emission (a drift flux) due to neutrino-nucleon scatterings is absent at the leading level 
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of 0(hbB/T), where \xb is the nucleon magneton, B is the magnetic field strength, and T is the 
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matter temperature at a neutrinosphere. This is because at this level the drift flux of the neutrinos 
are exactly canceled by that of the antineutrinos. Hence, the relevant asymmetry in the neutrino 
emission is suppressed by much smaller coefficient of 0(\ibB /m), where m is the nucleon mass; 
detailed form of the relevant drift flux is also derived from the scattering cross section, using a 
simple diffusion approximation. It appears that the asymmetric neutrino emission is too small to 
induce the observed pulsar kicks. However, we note the fact that the drift flux is proportional to 
the deviation of the neutrino distribution function from the value of thermal equilibrium at neu- 
trinosphere. Since the deviation can be large for non-electron neutrinos, it is expected that there 
occurs cancellation between the deviation and the small suppression factor of 0{^bB /m). Using 
a simple parameterization, we show that the drift flux due to neutrino-nucleon scatterings may be 
comparable to the leading term due to beta processes with nucleons, which has been estimated to 
give a relevant kick velocity when the magnetic field is sufficiently strong as 10 15 -10 16 G. 
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I. INTRODUCTION 



A core-collapse supernova explosion is one of the most spectacular events in astrophysics; 
99% of its gravitational binding energy is released as neutrinos, while only 1% as the ki- 
netic energy of a shock wave. Therefore, neutrinos play an essential role in supernova 
explosions, and their detection by ground-based large water Cherenkov detectors, such as 
Super-Kamiokande and Sudbury Neutrino Observatory, would provide valuable information 
on the nature of neutrinos as well as supernova physics (see Ref. and references therein). 

Because of their dominance in supernova energetics, the neutrinos may also give the 
solution to a long-standing astrophysical mystery, or "pulsar kicks." Recent analyses of 
pulsar proper motion give a mean birth velocity 200-500 km s , with possibly 

a significant population having velocity > 1000 km s _1 . These values are much greater 
than the velocities of their progenitors (~ 30 km s _1 ). A natural explanation for such high 
velocities is that supernova explosions are asymmetric, and provide kicks to nascent neutron 
stars. In this paper, we are concerned with models in which the pulsar kicks arise from 
magnetic field induced asymmetry in neutrino emissions from proto-neutron stars. We do 
not deal with another type of mechanism which relies on hydrodynamical instabilities in 



the collapsed stellar core 111 ], leading to asymmetric matter ejection and/or 



asymmetric neutrino emission. Concerning the latter mechanism, several numerical results 
gave rather negative results P. 11^. 

There are many past studies which have noted that parity violation in the weak interaction 
can lead to asymmetric neutrino emission in the strongly magnetized neutron star matter 
see, e.g., Refs. Q,Q,Q)- 

Arras and Lai 14] have shown that an asymmetric neutrino 
emission (a drift flux) that may give the pulsar kicks is mainly induced by charged-current 
interactions with nucleons, v e n — > e~p and u e p — *> e + n, in magnetic fields. On the other hand, 
the asymmetric emissions due to neutral-current interactions, vN — > uN, are found to be 
irrelevant because their contributions cancel between neutrinos and antineutrinos. Further 
their finding is that the drift flux exists only when the distribution deviates significantly 
from the value of thermal equilibrium. This result has invalidated the previous optimistic 
estimation that sufficient asymmetry can be obtained by multiple scatterings of neutrinos 
by nucleons, slightly polarized by the magnetic field of moderate strength (~ 10 12 G) jl^ . 
Their naive estimation has shown that the asymmetry parameter for the u e , v e flux would be 
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sufficient for the observed pulsar kick velocity, if the magnetic field is very strong, B ~ 10 15 - 
10 16 G. 

In the present study, we focus on the effect of weak magnetism and recoil corrections on 
neutrino-nucleon scatterings (particularly u^-N and v^ T -N) in supernova magnetic fields, 
which has not been considered in the previous publications such as Ref . Q| . We expect that 
including these two corrections changes the previous conclusions, which are listed above, as 
follows. The weak magnetism correction reduces antineutrino-nucleon cross sections com- 
pared to those of nentrino-nncleon scatterings Q. For the reason, the cancellation of 

n 

asymmetry between u and u, which were shown in Ref. |14j |. does not occur at 0{k/m) 
level, where k is the neutrino energy and m is the nucleon mass. Although the drift flux due 
to u-N scatterings is suppressed by a small factor of 0(k/m), this term may be as large as 
the term due to charged-current interactions between electron (anti) neutrinos and nucleons. 
This is because the drift flux is proportional to the deviations of distribution functions from 
the value of thermal equilibrium at a neutrinosphere Q] 1 For fi and r neutrinos, with 
which we are concerned in this study, the transport opacity is primarily resulted from u-N 
scattering, whereas energy exchange is due to elastic u-e~ scattering whose cross section 
is substantially smaller. As a consequence, the decoupling layer locates at much deeper 
region than the neutrinosphere and there should be large optical depth for the asymmetric 
flux to develop. On the other hand, for u e and u e , the difference from equilibrium at the 
neutrinosphere is considered to be very small owing to dominant u e n — > e~p and u e p — > e + n 
reactions, which act as very efficient energy-exchanging and thermalizing reactions. Thus, 
although the drift flux due to u^ T and u^ T is suppressed by a small coefficient of 0(k/m) 
compared with that due to u e and u e , we believe that the large deviations from equilibrium 
for non-electron neutrinos can counter the suppression factor. 

Motivated by the above reasoning, we first derive the cross section for u-N scatterings in 
magnetic field by including the weak magnetism and recoil corrections of scattered nucleons. 
In general, it is quite difficult to evaluate. However, in the limit that the nucleons are 
nondegenerate, we can present a simple form of the differential cross section, and give 
a diffusion equation which enables us to discuss qualitatively whether the effect we are 



1 In this paper, we use a term "neutrinosphere" as that above which neutrinos freely stream out, while 
by terms such as "decoupling surface," we intend the surface within which neutrinos are kept in thermal 
equilibrium. 
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tackling is important or not. We show, as a result of the diffusion equation and a simple 
parameterization, that the drift flux on which we concentrate can contribute at a considerable 
level for reasonable parameter choice. 

Throughout this pap er, we do not consider the effect of neutrino oscillation, although 
recent atmospheric |l7^ . solar jl^, and reactor experiments [19] have shown that the neu- 
trinos have the nonzero masses and mixings. This treatment is justified because the flavor 
mixing between electron and the other flavor neutrinos is strongly suppressed by matter 
effect. Roughly speaking, the efficient flavor conversion takes place at a so-called resonance 
region, where a condition n e ~ Am 2 /GFk is satisfied, however, at much deeper region 
n P p> Am 2 /Gpk matter induced neutrino oscillation is highly suppressed (see, e.g., Ref. 
|2(|). With the recently inferred parameters Am 2 by the oscillation experiments, the reso- 
nance region locates at rather outer envelope such as O+C or He layers; at a region around 
the neutrinosphere on which we are focusing in this study, matter potential strongly pre- 
vents electron (anti) neutrinos from mixing with the other flavor (anti) neutrinos. In other 
words, at the sufficiently deep region in stars, electron (anti) neutrinos propagate as the 
mass-eigenstates. For /i and r flavor (anti)neutrinos, they are found maximally mixing with 
each other, however, it does not affect the supernova dynamics such as pulsar kicks, because 
they can not be distinguished in supernovae. Thus, there is no need to worry about the 
effects of neutrino oscillation on supernova dynamics. Several years ago, Kusenko and Segre 
[2ll | proposed a very interesting mechanism that the pulsar kicks may be induced by neutrino 
oscillation. Unfortunately as we have already noted above, this mechanism does not work 
unless we introduce very heavy sterile neutrinos, because with the experimentally inferred 
parameters neutrino oscillation is strongly suppressed near the neutrinosphere. 

The remainder of this paper is organized as follows. In Sec. |Hj we derive a general formula 
for the matrix element and the cross section of the v-N scatterings by considering the weak- 
magnetism and recoil corrections. A more concrete form of the cross section is derived in 
Sec. based on the reasonable assumption that the nucleons are highly nondegenerate. 
In Sec. IIV[ using a simple diffusion approximation, the drift flux due to v-N scatterings is 
given in a very simple form, and we discuss the implications of the diffusion equation in Sec. 
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II. MATRIX ELEMENT AND DIFFERENTIAL CROSS SECTION 



We consider neutrino-nucleon scatterings via the neutral-current interaction in magnetic 
fields. The usual V — A current is modified when the weak magnetism correction is included 

as 



2m 



(1) 



where = — (k^)' is the momentum transfered to the nucleon. We label the neutrino 
momentum by k, and the nucleon momentum by p; for the final state momenta, the symbol 
' is attached. We summarize the values of couplings, cy, ca, and F 2 in Table |U 

TABLE I: Coupling constants. Here gA ~ 1.26, sin 2 9w ~ 0.23, [i p = 1.793, and \i n = —1.913. 



Reaction 



ca 



F-2 



Vp — > Vp 2 

vn — > vn 
v e n — > e~p 
v e p — > e + n 



2sin 2 6'iy ^ \{ftp — ftn) — 2 ftp sin 2 9\y 



2f ~\{ftp - ft n ) - 2ll n s\VL 2 9w 



OA 
9A 



ftp ftn 
ftp ftn 



Since the weak magnetism correction is at k/m level, we must also include the nucleon 
recoil correction of the same order. The initial and final states of nucleon spinor can then 
be expressed by 



u s N (p)u s N {p) 
u s n(p')u s n (p) 



1 - S7 5 7 3 ) ( 1 + 7° + — Y ) - *-7 5 (l + 7° 



P3 



m J m 
(1 - S ' 7 V) ( 1 + 7 ° + &>A - s'%\\ + 7 ° 



(2) 



m ) m 

where s, s' = ±1 specify the initial and final nucleon spins relative to the magnetic field axis 
(3- axis). 

With these corrections, the matrix element for u-N scattering is given, to the order of 
k/m, by 



\M ss >(n,n', p, P ^ 12 



M { 3{n, n') " + 6\M ss >{n, n', q)\' 2 + 5\M ss >(n, o', p)\\ (3) 



where Q = k/k and O' = k'/k'. Ai s J (f2, ft') is the leading term already given in Ref. 
Qby, 



m { 3(W) 



2 G 2 

= -if[<v + 3 4 + (4 - • n ' + 2 °a{cv + c A )(sQ + s'ft') ■ B 



+ 2c A (c v - c A )(sQf + s'Sl) ■ B + ss'{(c 2 v - c\){l + ft • SI') 

+ Ac 2 A n ■ BQ' ■ B}], (4) 

5\M. ss '(fl, ft', q)\ 2 is the correction term at 0(k/m) which depends on q(— p' — p = k — k') 
as well as f2 and fl f : 

5\M ss ,(n,n',q)\ 2 = ^[±2c A (c v + F 2 )q-(n-n')+2c A (c v + F 2 )(s-s r )q-B 

± 2c A (c v + F 2 )ss'q ■ B(n' — ft) B 
+ c A (c v + F 2 )(s' - s)(q ■ Sin' B + q ft ft ■ B) 
T c v {c V + F 2 )(s + s'){(fi' xOjxq}. B], (5) 

and S\M. ss '(ft : f2', P)\ 2 is also the correction term at 0(k/m) ) which depends on P(= p+p'): 

5\M s An, Of, P)\ 2 = ^[-{c 2 v + c 2 A + (c 2 v - c 2 A )ss'}P • (fi + fi') - 2c v c A (s + s')P • B 

2m 

- 2c 2 A ss'p ■ B(n + n')-B- c v c A (s + s')(p nn' b 

+ P ■ fl'ft ■ B) ± c A (s' - s){(fi' x n) x P} ■ B]. (6) 

Here, upper and lower signs represent the expressions for neutrinos and antineutrinos, re- 
spectively; this notation is kept throughout this paper. Note that we can explicitly check 
the time-reversal symmetry for the matrix element \A4 SS >(Q, f2',p,p')| 2 , by exchanging all 
the initial and final state quantities, i.e., fl <-> Q', P — > P,q — > — <j, and s <-> s'. 

The differential cross section (per unit volume) can be obtained from the matrix element 
through the expression, 

x \M,An,(i',p, P ')\\ (7) 

where //v = /n(E), f' N = /at (£")> /n(E) is the nucleon distribution function with energy 
given by 

exp[(£ - fj, N )/T\ + 1 

and is the nucleon chemical potential (excluding rest mass). In evaluating the momentum 
integral, we neglect the Landau levels of nucleons, and therefore the nucleon momentum is 
a well-defined quantity. This is justified because many Landau levels are occupied for the 
conditions in a protoneutron star, and the change in the available phase space due to the 
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Landau levels is negligible. By substituting Eqs. (HJ)-© into Eq. (JJJ), the u-N cross section 
in the magnetic fields can be calculated to the order ofk/m. Equation (JJJ) is written in more 
convenient form as 



dT 



dk'dQ' 



Ao(k, k', y!) + 6A±(k, k', fi')Q ■ B + 6A T (k, k', fj,')Q' ■ B 



+ A™ M (k, k , //) + 5B WM (k, k , fx')Q ■ B + 5C WM (k, k , ■ B, (9) 
where fi f = Q ■ Q' . A and 5A± are the terms which include the matrix element 



A (k,k',fx') 
5A±n ■ B + 5A T Cl' ■ B 



k 



12 



(2tt) 3 

k' 2 
(2^)3 



Y,\Mf)(n,n')\ 2 s ( qo ,q), 

s,s' 

J2\M { 3(n,n') 2 5s ss ,(qo,q), 



(10) 

(11) 



where S (q ,q) and 5S ss >(q ,q) can be derived from the "nuclear response function," 
SsAQo* Q) = S o(Qo, q) + SS ss >(q , q), defined by 

d 3 p d 3 p' 



S ss '(q ,q) 



(2nY5\p + k-p'~k')f N {l-f N ). 



(12) 



(2tt) 3 (2tt) 3 

So(qo,q) is the leading term when B = 0, and SS ss '(qo,q) is the correction arising from 
nonzero B, which can be simplified, for small B compared to temperature (/j,bB/T C 1), 
to be 



, , m 2 T 1 , ( 1 + e~ x ° 

bo{qo,q) = - — m 



8S ss i 



2ixq 1 — e 
m 2 T 



v l + e~ x °- z / 
5x 



with the parameters given by 

5b 

Z T 



Xq 



5x 



HbB 
2T 



2irq (e x ° + 1) (1 + e -3 *-*) ' 

(go - q 2 /2m) 2 _ /^v 
4T(q 2 /2m) T ' 
\ + 2mo\ s + A _ 2mgo s , 



(13) 
(14) 



5 



5 



(15) 



■ Q, 



The detailed derivations of these relations have been already discussed in Ref. [14], and we 
refer the reader to the literature. The terms including A^ M , 5B WM , and ^C 1 ^ 1 , which are 
smaller than A$ and 5v4± by factor k/m, are calculated by 



A WM + 5jB ^M fi . ^ + ^Wfj/ . ^ 



,12 



+ 



Y, [8\M.Ato, n', q) | 2 {S (g , g) + 6S ss r(q , q)} 



d 3 p d 3 p' 



(2n) 4 6\p + k- p' - k')f N (l - f' N )5\M ss ,(Sl, W, P) 



(16) 



(2tt) 3 (2tt) 3 

where in the second term, we cannot take the matrix element out of the integral because it 
depends on p and p f , unlike that in the first term which depends on q = k — k'. 

Equations (fTU|) . (JUJ), and (fTT)|) are quite difficult to evaluate at this stage. However, 
we can calculate it assuming that the nucleons are nondegenerate, as we show in the next 
section. In fact, this assumption appears to be appropriate, because the asymmetric drift 
flux can develop only when the neutrino distribution deviates from thermal equilibrium; this 
occurs in the regime where nucleons are nondegenerate. 



III. NONDEGENERATE NUCLEON LIMIT 



For nondegenerate nucleons, the characteristic neutrino energy transfer in each scattering 
is of order go ~ k(T /m) 1 / 2 <C k. The cross section sharply peaks around k' = k, and we 
can evaluate Aq, 5A±, A^ m , SB WM , and SC WM , in a series of the small parameter (T/m) 1 / 2 . 
We define the dimensionless quantities 

"4(1 - ^T 1 1/2 



m 



u 



k'-k 
~lk~' 



(17) 



where the range of u is from — 1/e to oo. Using the expansions of the momentum integrals 
given in Appendix [XJ we obtain the following expressions, 

(go - g 2 /2m) 2n 



So(qo,q) 



m 2 T ( 2vr 3 \ 2 
2ixq \m 3 T 3 J ^ 



4Tg 2 /(2m) 



A 
5A± 



a wm 



k' 2 



l-l(l + ^) eu + eu* + 0(e 2 ] 



,/2 



c 2 v + 3c 2 A + (c v - 4)/i'] So (g , q). 
mg ' 



(2k) 



-±G 2 F c/-f 



Cy ± 2c a- 



So(qo,q), 



5B 



■ WM 



k' 2 AG 2 F 
~ (2tt) 3 ^T 
x S (q ,q), 
k' 2 2G 2 F fi B B 
(2tt) 3 m f~ 



(4 + c\)^fq • (n + n') ± c A (c v + F 2 ) q ■ (n' - n; 



4c A F 2 ^k-2c A (2c v + F 2 )^q.n> 
q l q l 



(18) 
(19) 
(20) 

(21) 



mT f mq\ 
q 2 \ Tq 2 



± 1-^ -C V (cy + F 2 ) }q.n> 



S (q ,q), (22) 



(2vr) 3 m T 



-4c A F 2 ^fk'-2c A (2c v + F 2 )^f q .n 
q 2 q 2 



So(»,9), (23) 



q 2 \ Tq 2 

where all of them are expressed by u and e instead of go and q, through the relations such as 
qo = k — k' = —euk and q = k — k' = kfl — k(l + eu)Q'. In deriving the expressions above, 
we have used the equation e MAf / T = n(27r 3 /m 3 T 3 ) 1//2 that is valid when B — 0, to relate 
the nucleon chemical potential to its number density n (the corrections due to finite B 
are of order B 2 ). These expressions are valid under the conditions T <^ m, k (mT) 1 / 2 , 
HbB < T, and fc > k min = ii B B(m/T) 1 / 2 ~ IO^^IBuT- 1 / 2 MeV. All these conditions 
are satisfied for the case of interest in our study. Unlike the previous publication [14], we 
must keep in the bracket of So to 0(e A ) terms although we did not explicitly show them in 
Eq. (|18j). This is because they are necessary for the evaluation of the drift flux which is 
suppressed at Oikjm) level compared with the leading term. 

Thus, we have derived the useful form of the differential cross section of neutrino-nucleon 
scatterings, which enables us to discuss the parity violating effect of the neutral-current 
weak interaction related to pulsar kicks. At first sight, however, these equations are still 
too complicated to obtain an intuitive implication for the asymmetric neutrino emission in 
the supernova magnetic field. Therefore, in the remainder of this paper, we derive a simple 
diffusion equation with which we can see how large drift flux is obtained. More delicate 
treatment using numerical simulations is beyond the scope of this study; it is slated for the 
future work. 



IV. DIFFUSION EQUATION FOR A DRIFT FLUX 



The Boltzmann equation for neutrino transport is written in the form 



dfu(k) 
dt 



dfu{k) 
dt 



(24) 



where we include the collision term due to scattering alone; effects of absorption are not 
considered in this section. The collision term in Eq. can be written by the following 
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integral [14]: 



\dfvW 


= 1 dk ' 


[ dt 


sc 



poo p (TP 

Jo dk 'J m 'dk^V k ' )5f » + jD(A: ' k ' )5U 



(25) 



where 8f u and 8f' v denote the deviation of neutrino distribution function from the value of 
equilibrium, i.e., 8f v — f v — fjf 3 '. The coefficients of 5f u and 8f' v are 



C(k, k') 
D(k, k') 



-90 IT, 



1 - /f) + /, 



(0) 
v i 



-qo/T f(o)/ 

J v 



f, 



(0)1 



(26) 
(27) 



The nonlinear terms which are proportional to 8f v 8f' v were dropped since we consider the 
regime where the deviation from thermal equilibrium 5f v is relatively small (the regime 
where the diffusion approximation is valid). 

For the purpose of examining the macroscopic consequence of the asymmetric cross sec- 
tion, we expand the neutrino distribution function in spherical harmonics up to dipole order 
as 

5f v (k) =g(k) + 3£l-h(k), (28) 

where g(k) is the spherically symmetric deviation from the equilibrium, and h(k) represents 
the dipole deviation which leads to the flux. These quantities, g(k) and h(k) are related to 
the energy density and the flux through the relations 

k 2 dQ 



U v {k) 
F v {k) 



(2tt) 3 
k 2 dfl 



kf v = 



''" r f(?)(k)+g(k)]=U™(k) + ^g(k), 

J Z7T Z 



2tt 2 
k 3 



■h(k), 



(29) 
(30) 



(2tt) 3 j 2vr 2 

where the quantity with superscript FD represents the value in the case of thermal equilib 
rium. 

The first moment of Eq. yields 





\dfu(k)] 


/ 4tt 


dt 



V ± SVl + 5V^ M ± 8V^ M + 8V»> WM ± 8V% 



(31) 



where the subscript "0" of V and 5V represents that the terms are the same for neutri- 
nos and antineutrinos, while "±" that the sign of each term is different for neutrinos and 
antineutrinos. The term Vq contributes at the leading level, whereas the quantities with 
superscript "5," ¥1," and ll B, WM" are suppressed by factors of 0(fx B B/T), 0(T/m), 
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and 0(nB-B/m), respectively. As a result of lengthy calculations (details are given in Ref. 
each term is written as 



Vr 



svl 



3tt 

AG 2 F k 2 ^fi B B 
~T~ 



n(c v + 5c A )h, 



3tt 



-n- 



2/, 



(0) 



5V 



WM 



k 



2G 2 F k 2 

— n — 

37r m 



-3(4 + 34) 



9 + T 
1 



dg 



dk 
2/f 



B, 



(32) 
(33) 



+ 5(4 + 94) 



k 



+2(4 + 54) fc 



dk J 



h + 24 ( 1 - 2/<°> + ~— ) k^f + 2c\kT 



6T\ , dh 



k J dk 



: d 2 h 



S yWM 



5V 



B,WM 




lQG 2 F k 2 k 

n—c A [c v + F 2 )h, 

Sti m 

AGlk 2 a B B k 

-n— c A 



(34) 
(35) 



9tt 



T m 



c v (1 - 2/^) + (2c y - hF 2 y- \ k^- - c v kT 



-2c v + 5F 2 ) (l - 2/S°>) + 2cyA;^| 



dfc 2 



(36) 



5V 



B,WM 
± 



-n 



9tt 



Uk dk 2 ] 



69 (l - - 68A;^ 



,/ ^31(1-2^°))+^+ ^ 

<9A; 3 



2A; 69T 



k dk 



+2{3(l- 2/ f) + ^}^ + 4^T 



(37) 



Since the sign of the 0(^bB /T) term is opposite for neutrinos and antineutrinos, there 
occurs cancellation of the drift flux at this level. Hence, the leading drift flux due to u-N 
scatterings is suppressed by a very small coefficient of 0(fisB /m). 

When we discuss the terms which depend on quantities such as dh/dk, d 2 h/dk 2 , we use 
the lowest order expression for h, i.e., 



h(k) ~ -J_v(/W(A0+s 



3k 
2G 2 k 2 



' " n(cl + 5c 2 



3n 



(38) 
(39) 



In this and the following expressions for the diffusion equations, we drop the time derivative 
term; it corresponds to a rapid redistribution of matter temperature, whose time scale 
is of order the mean free path divided by the speed of light c, much smaller than neutrino 
diffusion time of the star. In addition, we assume that all the quantities except for h, such as 
g,fj)°\ are the same for neutrinos and antineutrinos. All these assumptions are considered 
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to be quite applicable, because the discrepancy caused by them is further suppressed to 
the extent in which we are not interested. With the above assumptions, we take the first 
moment of the Boltzmann equation (|24j) . and obtain the explicit expression for the neutrino 
plus antineutrino flux F v + F D : 

2 



F v {k) + F p (k) 



3k 



1 - 5 W ) VU u {k) 



+ 



eV>5U v (h) + e^k^8U v {k) + e^k 2 -^5U u (k) 



B, 



(40) 



where 6^ and 5^ are small numbers of 0(k/m), and e^ 2 \ and e^ 3 -* are those of 
0{hbB /m); they are given in concrete forms by 

k 1 



$<i) 



mcy + 5c\ 



dk 



6 (2) 



+2 (c 2 v + 5c A ) k 

4k c\ 
3m c v + 5c 2 A ' 
4k [i B B ca 
3m" T (c v + 5c\) 2 



-3(4 + 34) (l-2/f)-5(4 + 94)^ 



(41) 
(42) 



[2c v (c 2 v + 174) - F 2 (54 + 4) } (l - 2/J 



(o) 



-2cy(c 2 v + 5c 2 A )k d -^- 



(43) 



,(2) 



4/c \x B B 



(~'a 



3m T (4 + 54) 2 
-2c v (4 + { 
4\i B B c v c A 



. , 2T\ 2 / 34T\ 



2c v (4 + 54) /f - F 2 (54 + 4) - 



(44) 
(45) 



3m c v + 54 ' 

As expected and already pointed out by Arras and Lai |14(, the drift flux, which is along 
the direction of the magnetic field B, is absent at order 0(^bB /T); the leading term of the 
relevant drift flux is at 0(fisB /m) as shown in Eq. (|40p. 



V. DISCUSSION AND CONCLUSION 



Neutrinos carry away almost all the energy released by gravitational collapse of a super- 
nova core, ~ 3 x 10 53 erg, which is ~ 100 times the momentum associated with the spatial 
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motion of pulsars. Therefore, 1% anisotropy in the momentum distribution of the outgoing 
neutrinos would suffice to account for the observed kick velocities. From Eq. ([40)1 . the kick 
due to v-N scatterings is characterized by the fractional momentum asymmetry factor 



where R v (k) represents a radius of the neutrinosphere of energy k. To evaluate the value of 
a sc we must be careful, because the numerator of Eq. (}4l)j) contains 8U U which is difficult to 
estimate without numerical simulations. However, we believe, for non-electron neutrinos and 
antineutrinos this deviation 8U V can be large, which may cancel the small factors e^. This is 
because of the fact as follows. The electron neutrinos and antineutrinos are kept in thermal 
equilibrium slightly within the neutrino sphere by beta processes, v e n <-> e~p and u e p <-> e + n; 
leading very little difference from thermal equilibrium SU VejDe at the neutrinosphere. In 
contrast, for v^ T and v^ T the main source for the opacity is vN <-> vN reactions. However, 
since the nucleon mass is much larger than the relevant temperatures, energy exchange 
between neutrinos and nucleons is inefficient. Thus, the decoupling surface, within which 
neutrinos are kept in thermal equilibrium by the subdominant energy exchange processes 
such as iViV <-> NNuu, e + e~ <-> vi>, and ve <-> ue, locates at much deeper region than the 
neutrinosphere; it should lead rather large deviation SU u (k). 

While a detailed study using numerical calculations is not given in the present paper, we 
simply compare our new drift term in Eq. (f4*U|) with the previously derived drift term due 
to absorption Q|. According to Arras and Lai jlJ], the drift flux due to the absorption 
u e n — > e~p is described by 



a, 



sc 



J °°dkJRl(k)dn[F u (k) + F p (fc)] • B 
J™dkJRUk)dn[F u (k) + F p (k)]-r 

J °° dkRl(k) [e^6U v (k) + e^k-g- k 5U u (k) + e ^k%5U u (k) 




(46) 



JS° dkRl(k)i- Q 



F 



.drift 




e ahs 5U Ue (k)B, 



(47) 



■ (tot) *(abs) 

where = K n 



i ^ 

+ «0, K 



*(abs) 



(abs) 




[1 + e (Hu-k)/T^ with 




(48) 
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2 C A (C A + Cy) jJ Bn B 2 C A (CA ~ Cy) /I Bp B 



(49) 



c v + 3c 2 A T c v + 3c\ T 

In order to compare the difference between the drift flux due to u-N scatterings and that 
due to absorption, we define the quantity -R sc /abs by 



R 



sc / abs 



2JS°dk \e^5U^(k) + e^k§:6U^(k) + e^k 2 ^6U^(k) 



* (abs) 



(50) 



which is the ratio of the asymmetric fluxes due to the scattering and the absorption. The 
factor 2 in the numerator represents that the drift flux due to scatterings mainly comes from 
two-flavor (// and r) neutrinos, whereas that due to absorptions is only related to electron 
neutrinos. More precisely, R 2 (k) must be included in the integral, but for simplicity we 
neglect it. In order to evaluate -R sc /abs, we need a further assumption concerning the energy 
densities U Ve and U„ . Here we parameterize them using the "effective temperature" T e eff 



and Tf^ as follows: 



UuAk) 



k 3 



2vr 2 e k/Tf + 1 ■ 
k 3 1 

2ir 2 e k / T f + T 



(51) 
(52) 



where for further simplicity, T e and are assumed to be independent of k. As the other 
parameters at the neutrinosphere, we take T = 3 MeV, p = 10 12 g cm -3 , and Y e = 0.3, and 
the reactions (scatterings/absorptions) with neutrons alone are considered. 

Figure [U shows a contour map for the value of -R sc / a bs, which is projected against T e efT and 



T° . As discussed above, we believe that the value of T e is very close to the temperature 
at the neutrinosphere (3 MeV), on the other hand, we take Tff to the extremely high 
value, 8 MeV, though it may be unlikely to be realized. 2 This figure clearly shows that a 
significant part of the drift flux can be attributed to v^ T -N scatterings, for a considerable 
area in parameter space of (T e eff ,T^ ). In particular, they are comparable to the effect by 
the charged-current interactions, when T^ s is extremely close to the temperature at the 
neutrinosphere, and this situation may be realized actually. Thus, we cannot neglect the 
effect of u-N scatterings even though at the leading order it is canceled between neutrinos and 



2 Of course, taking too high value for the effective temperature violates the diffusion approximation, but 
we believe that it still gives some useful implications. 
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FIG. 1: Contour map for the value of -R sc / a b s projected against T* s and Tjf. 

antineutrinos; it may give considerable contribution to the asymmetric neutrino emission, in 
addition to the charged- current interactions. Whether we can actually obtain the observed 
pulsar kicks from this mechanism is beyond the scope of this study and it is a subject for 
the future numerical study. 

In this paper, we have derived the cross sections of u-N scatterings in strong magnetic 
fields, including weak-magnetism and recoil corrections. Since neutrino interactions are 
described by theory of the weak interaction which violate the parity, the asymmetric neutrino 
emission (drift flux) is induced, and it may give the origin of pulsar kicks. The drift flux 
due to u-N scatterings is suppressed at 0{\i B B /m) level, because the leading contribution is 
canceled between neutrinos and antineutrinos. However, since the drift flux is proportional to 
the deviation of neutrino energy density from the value of thermal equilibrium 5U U and it can 
be large for non-electron neutrinos, it can cancel the small suppression factor of 0(/ibB jm). 
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We have shown that the drift flux due to v-N scatterings may be comparable to the leading 
term due to beta processes with nucleons. It is expected to give the non-negligible effect for 
inducing relevant pulsar kicks. 



APPENDIX A: DERIVATION OF DIFFERENTIAL CROSS SECTION 



In this section, we derive differential cross section, Eqs. (|18j) - (|2Hjl . from more general 
form, Eqs. (JTUj). (JUJ), and (jTHjl. in the case of nondegenerate nucleons. 

We do not give the discussion for deriving the detailed form for nucleon response function 
S ss >, because it has been already discussed in Ref. 



and we basically follow their method 



in deriving the remaining part. Our interest here is concerned with the momentum integral 
which includes P, i.e., 



d 3 p d 3 p' 



(2ir) 4 6\p + k-p'-k')f N (l-f' N )P, 



(Al) 



(2tt) 3 (2tt) 3 

which appears in the second term in the right hand side of Eq. ()16p. For the term including 
p, we first use d 3 p' to integrate over 5 3 (p + q — p') and then integrate over the azimuthal 
angle for p, resulting in 



^- r dpp 3 f l dfi fi5(q + E- E')f N (E) [1 - f N (E') 
lixq Jo J-i 



(A2) 



2irq 

with fi = p ■ q/pq and E' = —fi B Bs' + (p + q) 2 /2m = —fi B Bs' + (p 2 + q 2 + 2pqfx)/2m. 
Care must be taken to correctly integrate over the energy-conservation delta function with 
coordinate //; the integral is only nonzero if the argument of the delta function is zero for 
jj, G (—1, 1). This condition, fi 2 < 1 is rewritten in the following expression for momentum 
p 

rg - q 2 /2m- fi B B(s - s')] 



2 \ 2 
P >Pm 



q/m 



Here we change variables from p to E = —s/j,bB +p 2 /2m in the remaining integral; it gives 



m q 

2nq 3 JE min 



dE 



<l«-^-^B{s-s<) 



f N (E)[l-f N (E + q )} 



mq 

q 2 



qo - 



2m 



f* B B(s - s') 



S SS '(qo,q), 



A12 



E n 



-fj, B Bs 



[g + q 2 /2m - fj, B B(s - s')] 
A(q 2 /2m) 



(A4) 
(A5) 
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Note that the last form of Eq. (jA4|) contains the integral which is also obtained from the 
integration of nucleon response function S ss >(qo,q). With the similar manner, the term of 
Eq. (|A1|) which includes p' is derived in the following form: 



where 



m q 

2irq 3 Je' . 

1 min 



dE' 



qo + 77 — hbB(s - s) 

2m 



f N (E' - g ) [1 - f N (E') 



mq 

q 2 



q 

q + ~ HbB{s - s) 

2m 



S S s<(qo,q), 



(A6) 



S ss i(q ,q) 



EL 



m 



dE'f N (E'-q )[l-f N (E% 



2nq Je' . 

j- mm 

R , , bo + g 2 /2m - v B B(s - s')} 2 



(A7) 
(A8) 



Aq 2 /2m 

From this point on, we show S ss > = S ss i for nondegenerate nucleons. We first define the 
dimensionless variables x' = (E' — Hn)/T and z = qo/T; the integral yields 

m 2 T r°° , , / 1 \ / 1 



S ss '(q ,q) 



2nq Jx' . 

1 min 

m 2 T 



dx' 



1/ U + e-*' 



2nq 1 — e 

Expanding x' in to linear order in B we find 



In 



1 + e ^in 

l + e 



+Z 



(A9) 



in 



6x' 



~\~ Soc . 
(g + q 2 /2m) 2 _ /^v 
4T(q 2 /2m) T 



-HbB 
2T 



1 + 



2mg 



_ 2mq , 



(A10) 



g / v q , 

Note that 5x' is identical to 5x [Eq. (|15))] which already appeared in deriving nucleon 
response function. For fe < 1, S ss i can be written as a sum of So, the zero field value, and 
SS ss i, the correction due to the magnetic field, i.e., 



S ss '(qo,q) 

So(qo,q) 

5S ss '(q ,q) 



So{q ,q) + 5S ss ,(q ,q), 



m 2 T 



2nq 1 — e 
m 2 T 



In 



f \ + e~ x o~ 
v 1 + e~ x ' 
5x 



2nq (i + 



e o 



e o 



(All) 
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In the limit of nondegenerate nucleons with /ijv/T <C — 1 and exp(/X7v/T) = (27r 3 /m 3 T 3 ) 1 / 2 n, 
we find 

1 , /l + e" :r o+^ 
In 



1 — e~ z \ 1 + e x o 
1 



1 + e x o~^) (l + e- x 'o^ 

(go - q 2 /2m) 2 



~ e x ° = exp 



(A12) 



T 4T(g 2 /2m) j ' 

This last form is the same as that of Eq. (B9) in Ref. jl4j ]. which is the S ss > limit for 
nondegenerate nucleons, guaranteeing the relation S ss r = S SS ' for nondegenerate nucleon 
limit. 

Hence, we obtain the expression for the momentum integral which includes P as follows: 



/ WfW?^ 4<y4(p + k - p ' - k ' )M1 " f ' N)P 

= —^-[qo-^BB(s-s')]S ss/ (q ,q). (A13) 



r 

The nuclear response function S , ss /(g , ?) can be easily expanded by e using the same method 
given in Ref. |l4| |: leading to Eq. (fTH|) . Thus, all the expressions of Eqs. (Jinj)-(j2HJ) are also 
obtained. 
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